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Abstract 

Polydisperse brushes obtained by reversible radical chain polymerization reaction 
onto a solid substrate with surface- attached initiators, are studied by means of an 
off-lattice Monte Carlo algorithm of living polymers (LP). Various properties of such 
brushes, like the average chain length and the conformational orientation of the 
polymers, or the force exerted by the brush on the opposite container wall, reveal 
power-law dependence on the relevant parameters. The observed molecular weight 
distribution (MWD) of the grafted LP decays much more slowly than the corre- 
sponding LP bulk system due to the gradient of the monomer density within the 
dense pseudo-brush which favors longer chains. Both MWD and the density profiles 
of grafted polymers and chain ends are well fitted by effective power laws whereby 
the different exponents turn out to be mutually self-consistent for a pseudo-brush in 
the strong-stretching regime. The specific values are, however, inconsistent with a 
standard self-consistent field theory of pseudo-brushes which predicts a much softer 
mushroom-like layer. 
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*email:milchev@ipchp. ipc.bas.bg 



1 



1. Introduction 

Due to their potential for practical application^, polymer brushes (layers of linear polymer 
chains end-grafted to a surface) have been the topics of vivid scientific interest. Most of 
the experimental and theoretical studies so far, have been devoted to monodisperse brushes 
of polymer chains in which the reaction of polymerization has been terminatedi^. The 
conventional way of forming such brushes with sufficiently high density, whereby pre-existing 
polymer chains with functionalized groups on the one end are adsorbed on a surface, turns 
out to be extremely slow and inefficientS. Thus, the attachment of polymer molecules, once 
the surface is significantly covered, is heavily suppressed by kinetic hindrance. 

An effective way to achieve sufficiently high grafting densities is to grow the layer in situ, 
monomer by monomer, from a functionalized seed carrying polymerization initiator^ Theo- 
retically, the growth of polymer chains from a surface was studied recently by Wittmer et al.i, 
combining elements of Diffusion-Limited Aggregation (DLA)I with the theory of polydis- 
perse strongly stretched polymer brushedH'i. Generalizing the "needle growth" proble m§, 
they view the formation of the brush as a particular case of Diffusion-Limited Aggregation 
Without Branching (DLAWB): An infinitesimally small incident flux of matter feeds the 
growth so that the polymer coils (in the good solvent regime0'0) can relax their structure 
in response to excluded-volume interactions with their neighbors while new free monomers 
are irreversibly attached at the active ends of the chains. This approach predicts strong 
differences between polymer molecules grown on surfaces and polymers formed in the bulk. 
The polydispersity of the former is much higher compared to that for the same reaction in 
the solution because in a brush the longer chains are more efficient in adding new monomers 
than the denser short chains which have to compete for fewer " random walkers" . Below we 
briefly sketch some of the main predictions of this approach (see Sec. p D| ). 

In the present work we try to extend the investigation of in situ grown polymer brushes, 
focusing on their properties under conditions of chemical equilibrium between the polymers 
and their respective monomers as sketched in Fig.^ To this end a fiat impenetrable surface 
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is densely covered with (fully activated) initiators. These sites serve as anchoring points for 
linear unbranched chains, so-called "Living Polymers" (LP)0, which can grow by reversible 
polymerization until an equilibrium with the ambient phase of free monomers is established. 
We suppose that both the scission energy J and the activation barrier B are independent of 
the monomer position (along the chain contour as well as in space) and density. Desorption 
events occur only at the active chain ends. In contrast to the systems of " Giant Micelles" 
(GM) considered elsewhere00, the grafted LP chains are not allowed to break and there 
are no freely floating LP in the bulk. The LP brush is grafted at the bottom (at 2; = 0) of a 
long container with an additional inert wall at the top {z = Lz), i.e. we consider a canonical 
ensemble and the total number of grafted and free monomers Mj = Ng + Nf is conserved. 

Although this system is qualitatively different from the case of irreversible DLAWB 
growthlll, some properties that we have investigated by means of an off-lattice Monte Carlo 
algorithm are found to match surprisingly well. Thus, the longer chains are again favored 
(compared to corresponding bulk systems) due to the gradient of the density profile of grafted 
monomers. Hence, the polydispersity again becomes very large and we observe (surprisingly) 
virtually the same power law behavior for the MWD c{N) as reported by Wittmer et al.i. 
This finding pertains to the established density profiles of aggregated monomers (f){z) and 
active chain ends p{z) as function of their distance z from the grafting surface as well. 

We will first sketch in Section ^ the elements of a theoretical description of a dense 
LP brush. The computational algorithm is then outlined in Section ^ The preparation of 
equilibrium configurations and the kinetics of layer growth until the onset of equilibrium is 
briefly examined in Section ^, and the main results of this paper on the static properties 
of LP brushes are presented in the subsequent Section |^. We examine the influence of 
total monomer concentration (pt and grafting density l/rf^ of the initiators on the density 
profiles, the average degree of polymerization (A^), and on the conformational properties of 
the polymer chains in the brush. In order to make connection with possible experiments 
using surface force machines, the force (F), exerted by the brush on the opposing wall of 
the container, is investigated. In Section |] we briefly discuss the computational results in 



view of the self-consistent field theory presented in Section and we conclude with a short 
Summary of the main results of this work. 

2. Analytical Considerations 
A. Mean-Field considerations 

In contrast to the DLAWB brush discussed by Wittmer at al.0, where the polymerization 
was assumed to be irreversible, we consider here a grafted polymer layer with annealed 
mass distribution in thermal equilibrium with ambient free monomers. Obviously, the task 
is to minimize the total free energy of the systemS which is in general a functional of the 
free monomer distribution u{z), of the density profile of the grafted monomers 0(2;) and of 
the MWD c{N). At low enough free monomer concentration we may write this total free 
energy as the sum of the free energy of the brush Fbrushy^iz) , c{N)] and the free energy of 
the free monomers Ffree[u{z), (piz)]. The latter contains the usual entropy of mixing term as 
well as an excess chemical potential due to the interaction with the dense LP layer formed 
around z = 0. From minimization of both terms, subject to the constraint of fixed total 
particle number Mt, the mean number of monomers aggregated in the layer (Mg) may in 
principle be calculated as a function of Mt and the interaction parameters. Being only 
interested in the functional form of the MWD and the layer profiles, we merely state that 
such an "equation of state" can be found and an arbitrary number of monomers Mg can be 
bound per unit surface (within some obvious limits). Hence, (ignoring fluctuations of Mg) 
we decouple the problem in two parts, one for the free monomers with constant chemical 
potential /i(Mf) = log(u(z)) + fiex{(p{z)) and one for the Mg monomers in the layer at same 
chemical potential. 

The simplified task is now to minimize the free energy of the brush subject to the 
constraint of Mg monomers per layer. It is natural to write this free energy in a Flory- 
Huggins like manner as 



4 



Ftrush[4>{z), C{N)] = c(iV) (log(c(iV)) + /iliV + F,ha^n[(^{z), c{N) , N]) . (l) 

Af=l 

Note that all densities are taken per unit surface and that the Boltzmann factor is set equal 
to one. The first term is the entropy of mixing for the grafted chains, the second entails 
the usual Lagrange multiplieiS for the conserved number of aggregated monomers Mg (see 
above). The last term is the free energy of a reference chain in the self-consistent density 
profile created by neighboring chains, which in general is a function of chain length and a 
functional of both (f){z) and c{N). The above mean- field equation assumes in particular that 
the length of a given reference chain is not correlated with the chain lengths of its neighbors. 
We note that all contributions to Fchain[4>{z) , c{N), N] which are linear in chain length (even 
if dependent on (j){z) and c{N)) can be incorporated in the Lagrange multiplier term ^iN 
and are hence irrelevant. 

Before we continue to tackle the full self-consistent field problem in See p E| within the 
strong-stretching assumption of polymer brush theory (Sec. |2 CD we need to recap some 
elements of the theory of polymers and polymer layers. In the reminder of this subsection 
we will consider some simple cases of the general problem formulated above where the 
free energy of the reference chain is a function of only. Hence, one readily obtains the 
equilibrium MWD from the functional derivation of Fbrush[c{N)] with regard to c{N) the 
equilibrium MWD 

ciN) CX expi-fliN - F,ha^n{N)). (2) 

Two particular simple cases where this applies are a forests of needles growing vertically 
on a flat substrate and Gaussian polymer chains (without excluded volume interaction) fixed 
at a surface. In both cases the free energy per chain is linear in N and the MWD is simply 
c{N) oc exp(-Ar/(A^)). 

We can now also address the problem of dilute non-overlapping LP (so-called mushrooms) 
in good solvent fixed on the surface. If the surface is penetrable the free energy of the 
reference chain is the same as for dilute LP in the bulk, i.e. Fchain{N) = (7 — 1) log(A^) where 
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7 = 1.165 > 1 due to long-range excluded volume correlations along the chain00. This gives 
rise to a Schultz-Flory distribution and was recently confirmed computationallylllE. The 
situation is slightly different for LP on an impenetrable surface which reduces the partition 
function and effectively repels the chains^. In this case one expects Fchain = —T\og{N) 
with the exponent r = 1 — 7^ where 7^ ~ 0.65 < 1. Hence, one expects to find the weakly 
singular MWD 

c{N) oc N~^exp{-ijN). (3) 
We will re-address this issueiil in Seel 



B. Weakly-stretched mushroom-like layer 

If the total amount of monomers Mg is now further increased and the mushrooms start to 
overlap, the description of the equilibrium LP layer becomes more complex. One expects 
the excluded volume interactions between the monomers to favor longer chains which can 
explore more dilute regions of the layer. It is hence conceivable that the MWD is not a pure 
exponential, but a power law (plus cut-off). 

There are two possible extreme scenarios or architectures which are relatively easy to 
compute. One is to suppose that the LP chains strongly overlap and form a strongly 
stretched LP brush. We will consider this in the next subsection. A second natural guess 
is to speculate that the excluded volume interaction is extremely strong (compared to the 
entropic terms in Fhrush) and that, therefore, the layer attempts to reduce chain overlap as 
far as possibleH. 

A geometrical construction for this limiting case is a self-similar mushroom-like structure 
with ^{z) oc z. Hence, every LP chain is essentially contained within its own excluded volume 
blob of size ^ containing g oc ^^^'^ monomers. Here, z/ ^ 3/5 denotes the Flory exponent^ in 
dimension D = 3. This forces the density of monomers to scale like (f){z) oc g{z)/C,{zY oc z~°' 
with a = 3 — 1/z/ ~ 4/3 and the density of ends like the blob density p{z) oc l/C,{zY oc z~^ 
with (3 = 3. From this the MWD is easily found using c{N)dN = p{z)dz. This yields 
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c{N) oc N~'^ with r = 1 + 2z/ ^ 11/5. Note that because a > 1 most of the layer mass is to 
be found at the lower cut-off should this scenario apply, i.e. if it indeed minimizes the total 
brush free energy. 

C. Strongly Stretched Alexander-De Gennes brush 

In this subsection we assume that the equilibrium LP layer is described by a strongly- 
stretched pseudo-brush of Alexander-De Gennes typei, i.e. a compact layer of blobs of size 
^{z) oc (f)(z)~^^^^'^^^\ Coarse-grained on distances larger than the blob size ^, the chain are 
described by "classical trajectories" 



where the monomer index n is counted from the grafted initiator [n = 1). Within this 
geometrical construction the chains are strongly stretched on distances larger than ^. Hence, 
the number of chains per unit surface at z is given by 



with e = 2z//(3z/ — 1) ^ 2/3, H denoting the upper edge of the pile of blobs while Nh = 
n{z = H) stands for the upper cut-off of the layer so that the ends of all chains passing 
through an arbitrary plane at z stay beyond that planeEl. For the last term in eq.@ we 
have used the fact that the number of chains equals the number of blobs at z (which itself 
scales like 1/^^) and the connection between cj) and ^ imposed by the fractal dilute structure 
within the blob. 

Anticipating our computational results on LP brushes (Sec.|^), we use power-law functions 
to express (f){z) oc p{z) oc z~^, c{N) oc N~'^ and z{n) oc n~'^-^, as in the case of the 
weak-stretching limit (Sec. f2 B| ). Plugging this into eqs. (0)-(|§), we obtain three independent 
equations for our four exponentsEl. We express them in terms of r, the exponent of the 




(4) 




(5) 



MWD: 
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(3zy-l)(r-l) 

1 + z/ - (1 - z/)r ^ ' 

= 1/(1 + a{l - z/)/(3z/ - 1)) = (1 + z/ - (1 - z/)r)/2z/ (7) 

(3z^ — llr — z/ + 1 
+ l + .-(l-.)r 

The second of these equations results from the "trajectory" eq.(H), the third from equating 
the first and last term in eq.(^. 

One possible solution of the above three relations is — and this might at first sight 
appear surprising — the set of exponents derived above for the weakly-stretched mushroom- 
like layer (with = u). This just means that there is a smooth matching between the 
regimes of strong and weak stretching — at least for the exponent Jll. 

Another possible set of exponents which is consistent with the three strong-stretching 
relations arises for a brush formed by irreversible diffusion-limited aggregation. 



D. Diffusion Limited Aggregation without Branching 

Here we briefly recall the results^ derived by combining a mean-field treatment of the dif- 
fusive growth with a scaling theory of the growing polymers. The theory considers the 
case of slow growth due to irreversible polymerization which is supplied by a constant in- 
finitesimal flux of incoming monomers. The continuum field description of the growth is 
based on the above strong-stretching equations and, in addition to this, to two standard 
equations for aggregation processes assuming a second-order process for the reaction and 
mass conservation in the "adiabatic approximation"!. It is possible to solve these non-linear 
partial differential equations by means of a power law scaling ansatzi. To be consistent, the 
Laplacian field of incoming random walkers, requires one additional equation between the 
exponents, (3 = 2, which fixes the above set of exponents in a unique way: r = 7/4, z/^ = 3/4, 
(5 = 2 and a = 2/3. Hence, the layer is much more brush-like with S,{z) oc ^/z than the 
fluffy mushroom-like layer considered above. These exponents, as well as the related scaling 
relations, have been reproduced by means of MC simulations showing the DLAWB to be 
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marginally mean-field like in 3D in contrast to the standard DLA proble: 



E. Dense LP brush 

Finally we consider the original problem of minimizing the free energy of a dense brush 
Fbrush[4'i^)^ '^i^)] limit of high coverage Mg. Following a standard SCF procedure!, 

we integrate out the (f){z) dependence of the free energy^ and assume the free energy for 
a given MWD c{N) to be well described by the Alexander-De Gennes blob construction. 
Hence, we neglect the logarithmic terms not taken into account in this steepest-decent 
approximation. 

The free energy of a reference chain still depends on the MWD and the mass distri- 
bution, obtained by functional derivation with respect to c{N), will be in general different 
from eq.(0). For the minimization it is more convenient to use the total excess energy 
in the layer due to the chain interaction Fev[c{N)] = J2n c{N)Fchain[c{N) , N] rather than 
Fchain[c{N), N]. We assume that Fev[c{N)] is well approximated by the total number of 
blobs within the layer 

F, „[c(Ar)] = / dz/^{z)^ oc / c/ri0(n)(2-+i)/(3-i) oc / dn { dn'c{n') (9) 

Jz Jn=0 Jn=0 \Jn'>n J 

where k = 1 + l/2v 11/6. For the last equation above we have used (t){n) = 
(^J^y^dn'cln')^' from the strong-stretching relation eq.(^). A careful computation 

of the functional derivation of the total free energy of the layer, eq.(|I|), generalizing eq.(0), 
yields 

= = log(c(Ar)) + ^^,N + dn dn'c{n')^ (10) 

where we have left out irrelevant constants and prefactors. The second term on the right 
side describes the (trivial) exponential cut-off. For thick layers (large Mg) we examine again 
the asymptotic power law behavior ai z <^ H. 

Anticipating the result of Fig.^ we may look for a solution of eq.(|TU|) with c{N) oc A^""^, 
that is, we require the integrand in eq.(|TO|) to be oc 1/n. This argument would not be 
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changed if additional logarithmic terms oc log(A^) are added by hand in eq.(|T0|)E3. It is 
simple to check that this yields r = — 1) + 1 = 1 + 2z/ = 11/5 which is the same 
exponent found previously in Sec. p B| for the mushroom-like layer with ^(z) oc z. Within the 
strong-stretching assumption eqs.(|) to (|) this yields the exponents i^_l = z/ = 3/5,/3 = 3 
and a = 2/3. 

This is obviously a peculiar result in that by using at various points the strong-stretching 
assumption we eventually obtain a fluffy layer. However, as was stressed at the end of See p C| , 
both physical pictures match smoothly and one might still regard this results as marginally 
self-consistent. 



3. MODEL AND METHOD 

Since we are not attempting to describe the properties of particular system (such as 
polystyrene, etc.) but rather wish to contribute to the general understanding of univer- 
sal properties of in situ grown LP brushes, we may use an algorithm in which each chain 
consists of coarse-grained monomers connected bonds representing at least, say, 3 — 6 chem- 
ical bonds along the backbone of a polymer chain. 

In the present investigation we have harnessed the coarse-grained bead-spring algorithm 
for polymer chainsEl already applied successfully for systems of GM in ref.0. Our description 
can therefore be brief. The main difference with regard to the previous study is that — 
following the model definition sketched in Fig.|l| — the monomers may attach or dissociate 
reversibly only from end monomers of grafted LP chains. 

Each bond is described by a shifted FENE potential where a bond of length / has a 
maximum at Imax = 1 

UfENe{1) = —K{lmax — hY ln 

where J corresponds to the constant scission energy introduced in the Introduction. Note 
that UpENEi} = ^o) = ~J and that Ufene near its minimum at Iq is harmonic, with K 



1 - 



- J 



(11) 
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being the spring constant, and the potential diverges to infinity both when / Imax and 
I — Imin = 2/0 — Imax- Following ref.lll we choose the parameters l^ax — = h — Imin = 0.3 
and K/T = 20, T being the absolute temperature. The units are such that the Boltzmann's 
constant ks = 1- 

The non-bonded interaction between effective monomers is described by a Morse-type 
potential, r being the distance between the beads 

UAiir) = exp[-2a(r - r^™)] - 2exp[-a(r - rmin)] (12) 

with parameters a = 24 and rmin = 0.8. Then the 6-temperature of the coil-globule transi- 
tion for our model is B ~ 0.62 so that at T = 1 we work under good solvent conditional. 

The model can be simulated fairly efficiently with a dynamic MC algorithm, as described 
previouslyEllli. The trial update involves choosing a monomeric unit at random and attempt- 
ing to displace it randomly by displacements Ax, Ay, Az chosen uniformly from the interval 
—0.5 < Ax,Ay,Az < 0.5. Moves are then accepted according to the Metropolis criterion 
and one Monte Carlo step (MCS) involves as many attempted moves as there are monomers 
in the system. In addition, during each MCS as many bonds as there are initiators in the 
system are chosen at random at the active ends of the grafted chains, and an attempt is 
made to break them according to the Metropolis algorithm. Attempts are also made to 
create new bonds between the end monomers and free monomers within the potential range 
of Um (i.e. a new bond with energy Ufene)^- 

In order to keep the system in equilibrium with the ambient phase of single free monomers 
and prevent the longer polymer chains from touching the top of the container we have used 
a rather low value of the bond energy J = 2. The lattice constant d of the square grid 
of activated initiators is taken as a rule as = 1 for the case of a dense brush and, as a 
special case of a loose mushroom-like layer, d = 4. Typically, boxes of size x x 
and periodic boundary conditions in x— and directions have been used in the simulations 
with Lx = Ly = 16, 32, 64 and 15 < < 256 (all lengths in units of Imax) for systems from 
Mt = 8192 up to 32768 monomers. 
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4. KINETICS OF BRUSH FORMATION 



Since the present work has been focused mainly on the equihbrium properties of in situ 
grown polymer brushes (Sec.|^), we briefly report here some salient features of the kinetics of 
growth by reversible polymerization at constant Mt. Some idea about the the conformation 
of the growing pseudo-brush may be gained from the series of snapshots of the polymer 
layer, Fig.Q, taken at successive times where the average chain lengthil 6 < (A^) < 20. 

At the start of the polymerization dimers (A^ = 2-chains) are predominantly created with 
a large number of unsaturated initiators present. Each newly attached monomer becomes 
then an active chain end itself (dark spheres in Fig.|^), and may attach other free monomers 
in its vicinity or detach itself from the chain. The snapshot on Fig.|^ corresponds to a 
configuration at equilibrium. Evidently, while the density of the brush immediately at the 
grafting surface is very high, the space above it is dominated by several very long chains 
which barely overlap, so-called mushrooms. While shorter chains appear to be somewhat 
stretched this is definitely not true for the long chains (see Sec. |5 C|) . 

Since the chemical reaction in our closed system proceeds at the expense of the available 
free monomers, it is clear that the observed kinetics of polymerization should differ from 
that studied at constant incident fluJ. Thus no clear cut scaling relationships have been 
observed with the time evolution of the quantities of interest, as for example, in the case of 
surface coverage with time shown in the inset of Fig.|^. We show there the fraction of grafted 
chains which are larger than N = 2 and the fraction of unsaturated initiators (i.e. chains 
of length A^ = 1) versus time. (The data has been gathered during runs of more than 10^ 
MCS and was averaged over 20 such runs each starting with different initial configuration.) 

We stress that the equilibration kinetics is rather slow compared to similar systems of 
GM in the bulkS which are allowed to break and where the dynamics is not limited to 
end monomer adsorption and desorption. This forced us to simulate relatively short chains 
compared to our previous studyllZl. 

The total number of monomers in the aggregate is initially diffusion-limited, i.e. {Mg) oc 
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V^, and then levels off in a non-trivial way (not shown). Similar graphs revealing no clear 
cut dynamic scaling have also been obtained for the time evolution of the mean chain length 
(A^) and the characteristic brush height {H) . In contrast, {H{t)) plotted against (A^) 
during the relaxation of the system to equilibrium shows a nearly perfect relationship as 

Eventually, in order to show that we have also reached equilibrium for the distributions, 
we show in Fig.^ the gradual formation of the typical density profile (j){z) with time t 
elapsed after the beginning of the polymerization. We see that gradually the power law 
profile (f){z) OC (plus exponential cut-off) is build up. As will be discussed below (see 
Sec. |5 D| ) the exponent is consistent with a ~ 2/3. 

5. EQUILIBRIUM PROPERTIES OF LP-BRUSHES 

Most of the equilibrium properties of the polymer brush are studied here with varying total 
monomer concentration (pt of the solution whereby the scission energy J = 2 (and other 
bonded and non-bonded energy parameters), the total number of monomers Mt in the box 
are kept constant as well as = Ly, whereas the box size in z-dimension is varied. We 
begin by analyzing the role of the box size on the fraction of grafted chains, the mean chain 
length (A^) and the force (F) measured on the opposite wall, before we proceed to analyze 
the MWD, the conformational properties and the density profiles in view of the theoretical 
questions formulated in Sec.||. 

A. Effects of box size variation 

We first examine the variation of the fraction of grafted LP chains (with N > 2) and the 
fraction of remaining unsaturated initiators with system density. As depicted in Fig.^ both 
are found to change linearly with (pf Note that these fractions do not add up to unity 
since there is always an important fraction of sites with only one monomer attached — such 
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= 2-chains are not counted as proper tethered polymers. 

The mean polymer chain length (A^) (i.e. the average degree of polymerization), the 
mean pressure (F) on the opposite wall at z = Lz and the mean squared bond length (/^) as 
a function of L^are depicted in Fig.^. For small boxes where < 80 the LP brush interacts 
directly with the opposite wall while for larger boxes this interaction is only mediated via 
the ambient free monomer pressure. 

In the small box limit (A^) grows with following a power-law relationship (A^) oc 4>t~^^'^ 
(full line in Fig.|^). This surprising result clearly underlines the difference of the present 
system with respect to living polymers in the bulk where (A^) oc Evidently this 

effect is due to the influence of the upper wall at distance on the polydisperse brush: 
as soon as Lz becomes sufficiently large the monomers begin to dissociate from the active 
chain ends and the degree of polymerization drops. 

The compression of the LP brush due to increasing confinement can also be clearly seen 
from the mean square of the bond length. It is evident that a steady shrinking of the elastic 
bonds with density takes place for < 80. 

In order to characterize the dependence of the pressure within the box on the total 
monomer concentration 0t, i.e. to obtain an equation of state, we have measured the force 
(F), exerted by the free monomers and the LP layer, on the top of the box. In order to 
measure the force exerted by the monomers on the upper wall directly, the wall has been 
supplied with a Morse potential, eq.(|12D, sX z = where we have kept the repulsive branch 
only. As expected, one observes a steadily decreasing {F) vs relationship with growing 
distance L^. For > 80 the total pressure on the opposite wall is mainly due to osmotic 
pressure of the free monomers in the solution, hence is proportional to the free monomer 
density u{z ^ H) oc l/L^. (See Fig.|TO[) In contrast, for < 16 the system attains its 
maximum density of (pt ~ 2.13 and the pressure rapidly grows as the hard core repulsion 
between the monomers comes into play. In between these extremes, however, there is an 
interval of heights Lz where the pressure is found to follow a power-law dependence on Lz 
with an exponent of ~ —3.5 (long dashed line). For these values of Lz the opposing wall is 
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in immediate contact with the polymer layer. 



B. Molecular Weight Distribution 

In Fig.^ we present the observed MWD c{N) (normalized to the mean number of chains of 
length A^) of grafted LP chains versus chain length A^. Systems with 16 x 16 x Lz containing 
8192 particles are considered here. Despite some scatter of data for the longer chains at 
higher concentrations, a power-law relationship c{N) oc A^""^ with r ^ 7/4 is evident. 
Hence, the MWD of grafted LP decays indeed much slower than the exponential behavior 
found in the corresponding LP bulk systems briefly mentioned in Sec. (|2 A|) . Evidently, this 
is due to the monomer density gradient within the dense pseudo-brush favoring the longer 
chains and is to be expected. Much more surprising is the value of the exponent found. It 
is clearly smaller than the exponent r = 11/5 expected for a mushroom-like layer and is 
(within numerical accuracy) identical to the prediction of the DLAWB problem described 
in Sec. p D| . This flnding clearly disagrees with the SCF approach presented in Sec. |2 E| where 



we estimated the free energy of the layer by counting the numbers of excluded volume blobs. 
The layer appears to be much more brush-like. This is the central result of the paper. 

We do not have yet any completely satisfactory explanation for the observed r. Pre- 
sumably this is caused by logarithmic corrections which lower the free energy with respect 



to our steepest-decent estimation of See p E| . The fact that the MWD is of power-law form 
shows that they scale like the number of blob term and cannot be neglected. 

Finite size effects in the present simulation are negligible. This is clearly shown in Fig.|^ 
where a much larger system of 32768 particles in a box of 32 x 32 x 256 with densely {d = 1) 
placed initiators on the surface reveals the same power exponent r = 7/4, as in Fig.^ 

If the distance between anchoring sites d is further increased and the chains become too 
small (compared to d) to overlap, we expect the distribution eq.(|^) for mushrooms grafted 
on an impenetrable surface with a weak power law c{N) oc A^~(i~'^«) with 7^ ~ 0.65 and a 
pronounced exponential tail for larger A^. This is clearly shown in Fig.|^ (spheres) for the 
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system 64 x 64 x 64 (i.e. of same total density (pt = 0.125) with large distance d = 4 between 
chain anchoring initiators. In the inset we test the power law decay anticipated for small 
due to the effective repulsion of the hard wall. Qualitatively, this is confirmed (dashed- 
dotted line). Because the chains within this regime are not very long one could not expect 
to obtain a better estimate. Essentially much larger (A^) und d are needed than we are 
at present able to simulate. Note that similar (finite-size related) problems have also been 
reported for the estimation of the 7 exponent from the MWD of reversible polymerization 
in the buHJi. 

C. Conformational Properties 

The observed power-law behavior of the MWD with r = 7/4 does not necessarily mean 
that we have a dense polydisperse brush in the strong-stretching limit. In order to check 
this assumption directly we investigate the mean squared end-to-end distance {Rl), the 
radius of gyration ^-Rg) and the mean squared vertical position of the end monomer (Zg). 
These quantities are plotted in Fig.^ as functions of the chain length A^. The straight lines 
with slope = 0.6 demonstrate that the few very long and, hence, non-overlapping chains 
(A^ > (A^)) behave, as one expects, like flexible coils in a good solvent. More interestingly, 
though, the shorter chains show a steeper increase with exponent z/_|_ ^ 0.75. (This is clearly 
seen from the end-to-end distance and the radius of gyration and to a lesser degree from 
(Zg).) The exponent = 3/4 is expected from eq.(|^) for a strongly stretched layer for a 
MWD described by exponent r = 7/4. Hence, the slopes of Fig.|^ and Fig.^ and the strong- 
stretching assumption, are consistent. Consequently, the new exponent coincides again with 
the DLAWB prediction, but not with our SCF approach. 

The stretching of the polymer chains in the LP brush can be also made visible by looking 
at the second Legendre polynomial P2{cos{9)) = [3(cos(^^)) — 1]/2 as a measure for the average 
orientation of the covalent bonds with respect to the 2;-axis, normal to the grafting surface. 
Here 6 measures the angle between the bond vector and the z-axis, normal to the grafting 
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surface. One would find P2 = —0.5 for bonds lying flat on the surface, P2 = 1 for bonds 
perpendicular to it, and P2 = if all the bonds are oriented completely at random in the 
3Z^-space. In Fig.^ we plot P2 versus the bond index n (between the monomer n and n + 1 
where the monomers are counted along the chain starting with the initiator n = 1) averaged 
over all chains (with n < N). One can see that the bonds between neighbor monomers 
stay predominantly vertical to the surface, and the smaller n, i.e. the closer the bond is to 
the surface, the more pronounced this orientation is. With increasing bond index n and, 
hence, decreasing monomer density (f){z{n)), the local orientations become more and more 
unrestricted by the presence of other chains. As shown in the inset, the data compares well 
with the power law fit P2{n) oc n"*^'^^. This underlines again the power law character of the 
LP brush. 



D. Density Profiles 

From the observed MWD and the strong-stretching behavior discussed above we may now 
anticipate power law decay for both the equilibrium density profile of end monomers p{z) oc 
and of grafted monomers (j){z) oc z~°'. Specifically, from t ^ 7/4 one expects /? ~ 2 and 
a ^ 2/3 where we have used eq.(P) and eq.(|D respectively. (Obviously, the same predictions 
could have been made starting from = 3/4.) 



This compares surprisingly well with the slopes shown in Fig.|10| and Fig.|Tl] for the dilute 
limit with Lz = 256. (Note that the a-value is slightly too large.) Hence, the observed set 
of four exponents r, z/_l, (3 and a is self-consistent. 

Here we have presented results for different concentrations with Mf = 8192 particles in 
a box of size 16 x 16 x Lz- All profiles are normalized to unity, e.g. /^f^ (f){z)dz = 1. We 
observe the same values of the exponents a and (3 when the total concentration is further 
decreased, i.e. Lz increased. Even at higher densities, up to (pt = 2.0 (at Lz = 16), both 
(f){z) and p{z) retain their power-law dependence on z. While the function p{z) remains 
largely unaffected, the effective exponent a decreases. This is expected qualitatively for a 
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brush compressed between two walls where at very high compression it becomes a melt of 
blobs of constant size ^. 

In Fig.|10| we have also included the density profile of the free monomers in equilib- 
rium, u{z). It becomes evidently constant for large z where the volume fraction of grafted 
monomers becomes negligible. Obviously, the chemical potential of the free monomers has to 
be constant with respect to position, i.e. /i = Hexi^z) + \og{u{z)) = const. Within mean-field 
one expects an excess chemical potential of fiex{z) oc (j){z) +u{z) which is ~ (f){z) for z <^ H. 
This is indeed born out as one verifies by plotting \og{u{z)) versus 0(2;) (not shown). This 
explains the slightly more complicated (not power law) form of the observed u{z). 



6. DISCUSSION 

To our knowledge, this is the first time that this problem (as defined in Fig.|T]) was considered 
by numerical as well as analytical means and, admittedly, our understanding of the LP brush 
is still incomplete. 

The effective exponents r = 7/4, z/_l = 3/4, /5 = 2 and a = 2/3 found in the simulations 
form a self- consistent set of exponents in agreement with the strong-stretching predictions 
eqs.(^||,^. Surprisingly, these values reveal a nearly perfect agreement (within inevitable 
statistical error) with the predicted exponent^l for the apparently different problem of the 
in situ formation of a grafted polymer layer through diffusive-limited aggregation with irre- 
versible polymerization (DLAWB). It is important to stress, however, that our values might 
still be somewhat off the DLAWB exponents due to the additional exponential cut-off (ef- 
fecting the long chains) and the too short mean chain lengths (A^) (hence, the too small 
number of decades) we have been able to simulate. 

Yet the question about the observed similarity between a pseudo-brush with an annealed 
MWD in thermal equilibrium and a steady-state DLAWB brush with a quenched MWD 
remains. In order to understand these results better we have attempted to describe the 
layer by a standard self-consistent field theory (SCF). As a steepest-decent estimation of the 
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free energy of our polydisperse LP brush we used the number of blobs contained in the layer 
ignoring all logarithmic corrections^. This yields a (marginally self-consistent) self-similar 
mushroom-like layer with ^{z) oc z. Hence, every chain should be essentially contained in its 
own blob. However, the resulting fluffy layer is inconsistent with the strongly overlapping 
pseudo-brushes with ^{z) oc ^/z we find in the computer experiment. In particular, the 
observed decay of the MWD with t ^ 7/4 is much slower than the SCF result r = 11/5. 
This is the surprising and central result of this paper. It clearly shows that the free energy 
of the brush in our approach has been overestimated and that the logarithmic corrections 
have to be taken properly into accountH. Unfortunately, these corrections are not well 
mastered even for monodisperse quenched polymer^l. At best one may take at present 
r as a phenomenological parameter characterizing the free energy of the layer. It would 
be interesting to estimate numerically the free energy of the layer by some direct method 
by varying the interaction parameters, although we have not attempted this yet. Thus 
further studies, numerical as well as analytical, would be required for a full insight into the 
equilibrium behavior of the LP brush. 



7. SUMMARY 

Polydisperse brushes of high grafting density, obtained by reversible radical chain polymer- 
ization reaction onto a solid substrate with surface-attached initiators, have been studied by 
means of an off-lattice Monte Carlo model of living polymers (LP). Basically, our simula- 
tional results for the layer in equilibrium indicate that a moderately stretched pseudo-brush 
is formed which is well described by diverse power law distributions plus an additional 
exponential cut-off at the edge of the layer: 

• The Molecular Weight Distribution (MWD) of the chains of contour length fits 
c{N) oc A^""^ with T ^ 7/4. Hence, the MWD of grafted LP decays much slower than 
corresponding LP bulk systems. This is due to the monomer density gradient within 
the dense pseudo-brush favoring the longer chains. 
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• If the chain overlap becomes small the MWD becomes essentially exponential with a 
weak singularity for small N with r — 1 — ^ 0.35 caused by the effective repulsion 
of the impenetrable wall. 

• The end-to-end distance and the radius of gyration show R oc N'^-^ with i>± pa 3/4 for 
smaller chains and i'± 0.6 for long enough, not overlapping mushrooms. 

• The density of active chain ends is well fitted by p{z) oc with j3 ^ 2 within a broad 
concentration regime. The density profile of aggregated monomers scales roughly like 
(j){z) oc z""" where the exponent becomes a ~ 2/3 in the dilute limit. 

Additionally, the influence of total monomer concentration and of grafting density of 
the initiators on the density profiles, on the average degree of polymerization, and on the 
conformational properties of the polymer chains in the brush have been studied too. We have 
also investigated the force, exerted by the brush on the opposing wall of the container. Very 
briefly, we have also examined the kinetics of brush growth until the onset of equilibrium. 
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Fig. 1. Sketch of our model of a hving polymer (LP) brush. Polymer chains grow reversibly from 
the (fully activated) starters (dark squares) fixed on the (impenetrable) bottom wall. They are in 
thermal equilibrium with a reservoir of free monomers (open spheres). We suppose that both the 
scission energy J and the activation barrier are independent of the monomer position (along the 
chain contour as well as in space) and density. Desorption events occur only at the active chain 
ends (dark spheres) and chains are not allowed break. Branching of chains is forbidden as well. 
The brush is put into a container with an opposite wall at z = L^- The total number of particles 
Mt is conserved. 
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Fig. 2. A snapshot of a brush configuration corresponding to mean chain length {N) 20. 
Neither the ambient free monomers nor the array of initiators are shown. Different intensity of 
grey color applies to different chains. Only the shorter chains do sufficiently overlap and, hence, 
are somewhat stretched by the density gradient. 
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Fig. 3. Activation of initiators at fixed total particle number Mt = 8192: Fraction of unsaturated 
initiators [N = 1: spheres) and of grafted cliains {N > 2: triangles) versus oc l/L^ (main figure). 
In the inset the same quantities are traced versus time elapsed after starting the polymerization 
reaction for the standard system with = 256 and (pt = 0.125. 
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Fig. 4. The density of grafted monomers (j){z,t) in a system of total concentration (pt = 0.125 
with Lz = 256 during relaxation to equilibrium. The resulting effective slope at late times is 
a 2/3. 
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Fig. 5. Effects of the box volume oc at Mt = 8192 and Lx = Ly = 16. The mean chain length 
(AT) (squares), the brush pressure on the opposite wall (F) (spheres) and the average squared bond 
length between neighboring monomers along the backbone of the polymer chains (diamonds) are 
plotted as function of L^. For large boxes > 80 the pressure is only due to the osmotic pressure 
of the free monomers and (F) decreases with the volume, i.e. with slope —1. The bond length 
becomes constant. For smaller boxes where the brush is directly interacting with the opposite wall 
we find (iV) oc L^-^ (full line) and (F) oc L^'^-^ (dashed line). 
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Fig. 6. MWD of a dense polydisperse brush vs N for different total concentrations oc 1/Lz with 
Mt = 8192 showing clearly a power law c{N) oc N^'^ with r « 7/4. This is the central result of 
this paper. 
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Fig. 7. MWD c{N) for two very large systems (both of total density (j)t = 0.125 containing 
Mt = 32768 particles) with high (d = 1 : triangles) and low (d = 4 : spheres) grafting density of 
the initiators. The first system confirms the power law behavior with r ~ 7/4 (dashed line) while 
the mushroom- like system shows for large essentially an exponential decay (full line). However, 
as can be seen from the inset where we have retraced the data in log-log coordinates, this MWD 
is not incompatible with the weak initial power law behavior (dashed-dotted line) expected for 
non-overlapping mushrooms fixed on an impenetrable wall. 
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Fig. 8. Radius of gyration (spheres), end-to-end distance (R^) ' (triangles) and mean 

squared end monomer position (z"^)^'^ (full diamonds) averaged over chains of given mass N. Note 
height {H) = 18.2 and mean chain length {N) = 33.9 of this configuration at (pt = 0.125. Dashed 
lines denote the slope of i'± = 0.6 corresponding to unstretched swollen coils in the good-solvent 
limit. The full lines indicate the slope u± = 0.75 consistent with a strongly stretched LP brush 
characterized by the exponent r = 7/4. 
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Fig. 9. Variation of the "order parameter" of bond orientation, P2{cos{9)), with bond index 
n. The data were obtained for a large box 32 x 32 x 256 containing 32768 monomers. The same 
data is depicted in the inset in logarithmic coordinates and compared with the power law fit 
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Fig. 10. Density profiles for the active end- monomers, p{z), in systems of Mj = 8192 particles 
within a 16 X 16 X box. The density of free monomers u{z) (stars) at Lz = 256 is included as 
well. The full line denotes the exponent f3 = 2 expected from the MWD and the strong-stretching 
assumption. 
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Fig. 11. 



Density profiles (piz) of aggregated monomers vs distance from the grafting surface. 



Same symbols and configurations as in Fig. IC. The full line denotes a slope of a = 2/3 consistent 
with eq.(P) and r = 7/4. 
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